Residence time and collision statistics for exponential flights: the rod problem 
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Many random transport phenomena, such as radiation propagation, chemical/biological species 
migration, or electron motion, can be described in terms of particles performing exponential flights. 
For such processes, we sketch a general approach (based on the Feynman-Kac formalism) that is 
amenable to explicit expressions for the moments of the number of collisions and the residence time 
that the walker spends in a given volume as a function of the particle equilibrium distribution. We 
then illustrate the proposed method in the case of the so-called rod problem (a Id system), and 
discuss the relevance of the obtained results in the context of Monte Carlo estimators. 



I. INTRODUCTION 

The so-called Pearson random walk describes the evo- 
lution of particles starting from a point-source and per- 
forming straight-line displacements until collision events, 
where cither the direction of propagation changes at 
random with probability p (scattering), or the trajec- 
tory is terminated (absorption) [TJ [2]. When the tra- 
versed medium is homogeneous, so that the scattering 
centers are uniform, the inter-collision distances (flights) 
are exponentially distributed. Exponential flights are key 
to understanding the dynamics of many transport pro- 
cesses, encompassing areas as diverse as radiation trans- 
fer, electron motion in semiconductors, gas dynamics, 
and search strategies [3HS]- In most such applications, 
one is typically interested in assessing the particle den- 
sity ^y (or some functional defined on $y) averaged over 
a d— dimensional volume V in the phase space. In Reac- 
tor Physics, for instance, ^y might represent the number 
of particles escaping from radiation shielding 4J . 

The particle density of exponential flights, in turn, is 
intimately connected to the statistical properties of the 
collisions ny falling in the region V , and the residence 
time tv spent within V [SHU]- Even under simplifying 
hypotheses, namely, that scattering and absorption prob- 
abilities do not depend on particle energy, so that we 
can safely define an average speed v (one-speed approx- 
imation), and that scattering is isotropic, the interplay 
between '^v, ny and ty turns out to be a deceivingly 
simple problem, and has attracted a renovated interest 
in recent years [SHE]- For instance, when the typical 
size R of the volume V is much larger than the mean free 
path At between collisions, namely, R^ Xt (the so-called 
diffusion limit), the normalized distribution of collision 
number V{ny) and the normalized probability density of 
residence times Q{ty) converge to each other (TUI, as il- 
lustrated in Fig. IT] (left). This is in general not true when 
R is comparable to At, i.e., when finite speed effects and 
boundaries come into play, and particles spend only a 
limited number of collisions in V before wandering away. 
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as shown in Fig. IT] (right) . 

Aside from its theoretical interest for understanding 
the dynamics of exponential flights in bounded geome- 
tries, the study of V{ny) and Q{ty) is also motivated 
by their prominence in Monte Carlo methods. In view 
of the intrinsic stochastic nature of exponential flights, 
one is naturally led to resort to Monte Carlo simulation, 
which can guide the development of analytical solutions, 
and, in most realistic applications, provide the answers 
that are not accessible by analysis alone [HI [17] . In plain 
Monte Carlo methods, the volume-averaged particle den- 
sity '^y for one-speed transport can be estimated by sim- 
ulating particle trajectories in the phase space and either 
counting the collisions ny in V ^ or measuring the length 
y of particle tracks within V . In the former case, we 
V~^\tY^i^yl = nyXt/V, whereas in the 
= y'^Y.^ev^^ = ^v/V [H]. As speed is 
assumedly constant, we can equivalently compute vp^"^'^' 
by measuring the residence time ty = iy jv that parti- 
cles spend in V . It is well known that the two estimators 
described above are unbiased with respect to ^y^ which 



have *^°" = 
latter ^^{^""^ 



amounts to saying that '^\ 



*v and * 



track 



*l 



where • denotes averaging with respect to particle trajec- 
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FIG. 1: The normalized distributions of collision number 
'P{nv) (dots) and residence times Q{tv) (solid line) in a vol- 
ume V for Id exponential flights. Left. _R = 10 and At = 1. 
Right. R=l and At = 1. 
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FIG. 3: Second moment of first-passage time (dark surface) 
and collision number (light surface) for leakage boundary con- 
ditions and p = 1. 



FIG. 2: Mean first-passage time and collision number for leak- 
age boundary conditions and p = 1. The two surfaces coin- 
cide. 



tories, and the limit is attained for an infinite number of 
realizations |17j . This in particular implies that the two 
estimators are related by fiv = W JTt, where Tt = Xt/v 
is the average flight time. Such a property non-trivially 
holds for any kind of boundaries imposed on V and stems 
from the memoryless (Markovian) nature of the under- 
lying exponential flight process [121 [H] . Fig. [T] suggests 
that tv and riy, while preserving the same average, will 
generally have different higher order moments, and in 
particular different variances. Hence, there might be an 
advantage in using either estimator for determining the 
desired particle density ^y. 

In the following, we address the issue of characteriz- 
ing the distribution of collisions ny and residence times 
ty in a volume V for exponential flights. This paper 
is structured as follows. In Sec. |llj we will first recall 
some preliminary background, and sketch a general ap- 
proach for the moments of the distributions, based on 
the Feynman-Kac formalism. Then, in Sec. |III| we will 
exemplify the proposed methodology by explicitly evalu- 
ating those moments for a Id domain, the so-called rod 
problem. Perspectives are finally discussed in Sec. |IV| 



II. METHODOLOGY 

The trajectory zj — {vt,ujt\ of exponential flights in 
the phase space is defined by the stochastic evolution 
of position and direction, starting from the point-source 
Zq = {ro,wo} at time t — 0. Due to the exponential 
nature of the displacement lengths, the stochastic pro- 
cess Zt is Markovian. In other words, knowledge of the 
pair position-direction at a given time enables to deter- 
mine the system evolution J37! . For the sake of simplic- 
ity, we assume that scattering is isotropic. The propa- 



gator ^(r, w,i|ro, wq) defines the probability density for 
the walker being at a point {r, w} in the phase space, at 
a time t, having started from the initial condition. The 
propagator of exponential fiights satisfies a probability 
balance, the forward Chapman-Kolmogorov equation 



d_ 
di 



*(r,w,i|ro,wo) =/:*(r,w,i|ro,wo), (1) 



where £ is the forward transport operator 
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(2) 



Here we have set v = cjw, t^ = As/w, As being the scat- 
tering mean free path, and the integral over directions is 
normalized to VL^i = 2tt'^/'^ /T{d/2), i.e., the surface of the 
unit sphere. We introduce then the collision density 



1 [+°° 
*(i",w|ro,a;o) = — / *(r, w,<|ro, wo)di, 
n Jo 



(3) 



which intuitively represents the equilibrium distribution 
of the particle ensemble. Remark that the propagator 
^(r, oj, t|ro, wq) depends on the boundary conditions im- 
posed on dV. The absence of boundary conditions corre- 
sponds to defining a fictitious ('transparent') volume V, 
where particles can indefinitely cross dV back an forth. 
On the contrary, the use of leakage boundary conditions 
leads to the formulation of first-passage problems p!5H2T] , 
where the walker is lost upon crossing dV. 

The collision number ny and the sojourn time ty of 
the walker inside V depend on the realizations of the 
trajectories Zj, and as such are random quantities, whose 
behavior can be fully characterized in terms of their re- 
spective moments 
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)(ro,a;o)= ^ nyV{ny\ro,ujo) 



nv— 1 

(i^)(ro,a;o)= / t'^Q{ty\ro,u;o)dty. (4) 
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FIG. 4: Mean residence time and collision number for trans- 
parent boundaries and p = 0.5. The two surfaces coincide. 



Remark that 7^(nv'|ro,a;o) and Q(iy|ro,a;o) depend on 
the initial conditions. Knowledge of all moments suffices 
to describe the associated distributions. In the following, 
we derive explicit expressions that allow evaluating the 
moments (n™)(ro,a;o) and {ty){ro,ujo) in terms of the 
equilibrium distribution $(r, aj|ro, wo)- 



A. Residence times 

In a series of seminal works based on Feynman path- 
integral formalism, Kac P^^ES] has worked out a general 
method for deriving the residence time distribution when 
the underlying stochastic process is a Brownian motion 
Wt, and later showed that his results hold more gener- 
ally for Markov processes [HS]. For a review (focused on 
Brownian motion), see, e.g., [27]. When a trajectory Zj 
is observed up to a time i, the associated residence time 
wit) < i in y is formally 



tv{t) = / x[^it')]dt', 



(5) 



x[z] being the marker function of the domain V , which 
is equal to 1 when z ^ V, and vanishes elsewhere. When 
V has leakage boundary conditions, tv{t) for an infi- 
nite observation time corresponds to the first-passage 
time to the boundary dV. More generally, the defi- 
nition in Eq. ([5| allows for multiple exits and re-entry 
crossings of dV [S] . The key ingredient of Kac approach 
is the stochastic integral F(t,s\zo) = (e"**^^*-*), where 
the expectation is taken with respect to the propagator 
^(z,i|zo), i.e., the probability density of performing a 
trajectory from zq at i' = to z at time t' = t, namely. 



Stv{t)\ _ 



)=yvI/(z,i|zo)e-^*^'(*)dz. 



(6) 
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FIG. 5: Second moment of residence time (dark surface) and 
collision number (lights surface) for transparent boundaries 
and p = 0.9. 



The existence and well-poscdncss of Eq. (l6| is discussed 
in, e.g., [22l |25]. By slightly adapting the treatment 
in [5S], it can be shown that F(t, s|zo) satisfies the equa- 
tion 



d_ 
di 



F{t, s\zo) = C*F{t, s\zo) - sx[zo]F{t, s\zo) 



where C* is the backward transport operator 



£* = Vo • Vr„ 



dujo — 



1 



(7) 



(8) 



Kac has shown that F(t,s\zo) can be interpreted as 
the Laplace transform (the transformed variable being 
s) of Q(ty|zo). The standard approach would therefore 
imply first solving Eq. ([7| for F(i, s|zo), and obtaining 
then Q(tv \zq) by performing an inverse Laplace trans- 
form. Eqs. (Tj and m are known as the Feynman-Kac 
formulae [57 . Once F{t,s\zQ) is known, the moments of 
residence time can be obtained from 



(f^')(zo,t) = (-l)' 



d" 



ds 



;F{t,s\za)\s=o. 



(9) 



Eqs. ([7| and ([9| yield the recursion property 
d 



dt 



(0(zo,t)=r(C)(zo,0 + "^x[zo](ty"')(zo,i), 



(10) 

with the conditions (i™)(zo,0) = and {ty){zo,t) = 1. 
In most applications, the observation time is assumed to 
be infinite, i.e., t — >■ -foo, which leads to the simplified 
equation 



r(C)(zo) = -mx[zo](t^'-^)(zo), 



(11) 



where tv{zo) = lim(_j.+oo ty(zo,i). Eq. ( [TTj ) is proposed 
in j8j to generalize a result by [11] and derive an elegant 
recursion formula for the moments of residence (and first- 
passage) times of exponential fights averaged over initial 
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conditions. In the context of Brownian motion, the rele- 




FIG. 6: Second moment of residence time (dark surface) and 
collision number (ligiit surface) for transparent boundaries 
and p — 0.5. 



vance of Eq. (11 1 is discussed at length in p8) 



When one is interested only in the moments of the 
distribution, and the observation time is infinite, the 
Feynman-Kac formalism may be rather cumbersome 



(Eq. (11) would still require inversing the backward op- 



erator C* ) , and can be altogether avoided by resorting to 
the so-called Kac moment formula [IS]. This approach 
has been successfully applied to the study of the resi- 
dence time of Brownian particles in |29j . For a review, 
see, e.g., |30| . The m-th moment of the residence time is 
obtained from 



(t^')(zo)=m! /dz. 



//>+00 I't2 

dzjn dtjn--- (iii*(z„,i,„ - t„_i|z,„_i) * ... * *(zi,ii|zo), 

Jo Jo 



where the convolution products read 



^(z.,r, t.,r - idzO * *(z. t. - .._,|z._,) = / ^z,*(z.,„ ..,, - t.|zO*(z., .. - t._r |z._r) 



(12) 



(13) 



Finally, by interchanging the order of integration in time, 
and extending the integration limit to infinity for each 
convolution product [H], we have the formula for the 
moments of the residence time, namely. 



(i")(zo) 



to! / dz„j*(z,„|z,„_i) * ... * *(zi|zo), (14) 



where 



1 f+°° 
*(z,+i|Zi) = - / *(z,+ i 

n Jo 



t\z.,)dt (15) 



is the collision density. Eq. ( 14 ) thus allows expressing 



the moments of the residence time as a function of the 
particle equilibrium distribution. 

B. Collision number 

In a previous work [TIP , we have explicitly derived the 
moments of the collision number ny for a broad class of 



renewal processes, when the point-source emits isotrop- 
ically. For exponential fiights, it is sufficient to remark 
that the process r„, i.e., the direction- averaged position 
of the walker, is Markovian at each collision event. The 
probability of performing ny collisions in the volume V 
is related to the propagator by 

V{nv\ro)^ /dr*(r,ny|ro)- /"dr*(r, ny + l|ro)(16) 



We introduce the direction-averaged collision density 



1'(r|ro) = ^*(r,n|ro). 



(17) 



n=l 



The derivation of the moments (ny)(ro) closely follows 
that of (ty)(zo) [TUJ. Here we just recall that the mo- 
ments of nv are given by 



-. m 

.')(ro) = - Vsm,fe/Cfc(ro), 



(18) 



fe=i 
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FIG. 7: Second moment of residence time (dark surface) and 
collision number (light surface) for transparent boundaries 
and p — 0.1. 



where the coefficients 



■ j=0 



(19) 



are the Stirling numbers of second kind [3T] , and 



Cfc(ro) = fc! / dvk... f dri*(rfe|rfc_i)...*(ri|ro) (20) 
Jv Jv 



are defined as fc-fold convolutions of the collision density 
*(r|ro) with itself [ID]. 

Now, the moments {ny){ro,uJo) for a directed source 
(5 (r — Tq) S {uj — ujq) can be evaluated as follows. First, 
we compute the density 7r(r'|ro,wo) of the walkers en- 
tering their first collision at r'. Each first-collision point 
will re-emit isotropically after the collision, i.e., the dis- 
tribution of the outgoing w' is uniform. Then, the mo- 
ments {ny){ro,ujo) are obtained by convoluting Eq. (18) 



for an isotropic source at r' with the first-collision source 
p7r(r'|ro, Wo). In terms of collision number probabilities, 
we have 



V{nv\ro,uJo)^P / dr'x[r']V{nv - l|rXr'|ro, ^o) +P / dr'x[r']V{nvW)n{r'\ro,oJo), 



(21) 



where xW] vanishes for r' ^ V and is equal to one else- where. This leads to 



){ro,ujo)=pYl (7)/^'"'^[^']<"^)('"')^(^>0'^")+^/^^'<O(^X''>0'^")- (22) 



r 



For the average collision number, from Eq. ( 22 ) we have 
in particular 



7r(r'|ro,^o)- 
(23) 



{nlr){ro,u}o) -= / dr' xW] +P dr*(r|r') 
By remarking that 

x(r)7r(r|ro,a;o)+p / dr'^(r|r')7r(r'|ro, wo) = ^{r\ra,ujo), 

(24) 
we therefore get 

(ni.)(ro,a;o) = / (ir^(r|ro, wq)- (25) 



Observe that for the mean residence time we have 



n 



dr / dw*(r,w|ro,a;o) 



(26) 



Since ^(r|ro,wo) = / da;^(r,a;|ro, wq), it finally fol- 
lows the (rescaled) equality between the average resi- 
dence time and the average collision number, namely, 

(n^)(ro,a;o) = (i^)(ro, wo)/Tt. 



III. THE ROD PROBLEM 

The approach presented in the previous Section al- 
lows explicitly evaluating the moments (n™)(ro,wo) 



<ny^>-<tv^> 




FIG. 8: The quantity (n\,){vo) - {tv)(xo) when xq 
function of 7? and the scattering probability p. 



= 0, as a 



and (iy)(ro,aJo)- When the equihbrium distribution is 
known, this amounts to solving the convolution integrals 
in Eqs. (22) and ( [l4| , respectively. However, analytical 
expressions for \I'(r7aj|ro,a;o) or ^(r|ro) (subject to the 
appropriate boundary conditions) are known only in a 
few cases O [TU], so that one must generally resort to 
numerical integration. A well-known and long-studied 
example where calculations can be carried out analyti- 



cally is the so-called rod model, where particles can move 
along a straight line [IJ [21 [32] ■ This corresponds to ex- 
ponential flights in Id, with only forward and backward 
direction allowed. Though the rod model is somewhat 
inadequate to address realistic radiation transport phe- 
nomena, we shall discuss it here for two main reasons. 
First, it allows illustrating the application of the above 
formulas for (n™)(ro,a;o) and (i™)(ro,a;o), and provides 
some hints on their use as Monte Carlo estimators. Sec- 
ond, the rod model, despite being admittedly oversimpli- 
fied, is nonetheless widely used in biology (often called 
velocity jump process), gas dynamics (Lorentz gas), fi- 
nance and neutronics, as it captures the essential features 
of the corresponding physical system [32H36] . 



We define ujf and LOb the forward and backward direc- 
tions, respectively. Similarly, we denote by Sf and Si, 
the forward and backward components of the source, lo- 
cated at a;o. Furthermore, we denote by x the abscissa 
of the rod, positive when oriented as w/. We set the 
mean free path Aj = 1, and we take v = 1. Scattering 
is isotropic. The volume V is assumed to be the interval 
[— i?, i?]. With this choice of parameters and notations, 
Eq. (IT]) reduces to the following set of stationary first- 
order differential equations 



d_ 
dx 



+ l\ '^{:x,ujf\xQ,Uf) = - [^{x,ujf\xQ,uif) -f ^(a;,a;t|a;o,w/)] + Sf 

d \ p 

— + 1 I 'i'{x,ujb\xo,ujf) = - [•^{x,ujb\xo,ujf) + 'i'{x,ujf\xo,ujf)] 



r 



(27) 



when the source is Sf, and 



— + 1 j 'i'{x,ujf\xo,ujb) = 2 [*(^''^/|2^o,Wb) + -^{x,uJt\xo,ujb)] 
— + 1 j 'i'{x,ujb\xo,'^b) = 2 [^{x,'^b\xo,uJb) + 'S{x,ujf\xo,ujb)] + Sb 



(28) 



r 



when the source is Sb- 

Two relevant examples will be considered here: i) 
leakage boundary conditions (a first-passage problem) 
without absorption, and ii) transparent boundaries with 
absorption. In the former case, leakages at x = ±i? 
impose "i'{—R,u!f\xo,ujf) — 0, ^(i?, ajb|a;o, w/) = 0, 
^(i?, Wfcjico, Wfe) — 0, and ^(— i?, cj/|xo,a;f,) = 0, which 
corresponds to an homogeneous medium surrounded by 
vacuum. The source is therefore ccq G V ■ In the latter, 
boundary conditions are imposed at infinity, which corre- 



sponds to an infinite homogeneous medium, the bound- 
aries of Y being transparent and not affecting particle 
trajectories. One-dimensional exponential flights are re- 
current walks (i.e., they almost surely re- visit their initial 
position) [HI [To] , so that it is necessary to impose leak- 
ages and/or set p < 1 in order to prevent (ny)(xo,a;o) 
and (iy)(xo,a;o) from diverging. 

For the case of purely scattering media, i.e., p = 1, 



and leakage boundaries, the rod problem equations (27) 



and (28) are straightforwardly solved by direct Integra- 



tion, and give rise to first-order discontinuous polyno- 
mials, the discontinuity being located at xq, i.e., at 
the source. Once the four solutions ^(x,aj^|a;o, wj), 
^{x,uJb\xo,ujf), '^{x,uJb\xo,i^b), and *(x, w/|xo, Wb) have 
been obtained, the moments (riy)(xo) and (i™)(xo) 
are computed by performing the convolution integrals 
in Eqs. (22 1 and (14), respectively. Remark that the 
isotropic source corresponds to assuming Sf = Sb and 



integrating with respect to the initial direction. For the 
mean first-passage time, we have 



(iy>(xo) = 



R^ + 2R^ 



with {nl,){xo) = {tv)i^o) (« 
the second moment, we have 



(29) 



1, so that Tt = 1). For 



("y)(xo) 



12R + 30i?^ + 20i?^ + 5i?* - 6x^ - URxf, - QR^xf, 

12 



(30) 



and 



(4)(xo) = 



12i?2 + 20i?3 + Si?-* 



12x1 



12Rxl - &R^xl + X 



12 



(31) 



The terms in these formulas look inhomogeneous (this 
is due to setting At = 1), but expressions are indeed 
dimensionless. The surfaces are discontinuous, since 
l^ol < R- Observe that when R is large we have 

(ty)(xo) ~ (ny)(xo). When R -^ +cxd, the moments (ny)(xo) > (t^)(xo) 
diverge, as expected from \d exponential flights being 
recurrent walks. Observe that the first and second mo- 
ment of the first-passage time satisfy the recursion prop- 
erty derived in [5], namely 



{(i™-i)(xo)}^ = 



{(iy)(xo)}5 



goes along with the collision and track length Monte 
Carlo estimators being unbiased with respect to each 
other. The second moments (ny)(xo) and (ty)(xo) are 
displayed in Fig. [3| It is immediately apparent that 
'actually, equality is attained only 
for i? ^ 1): this means that for this example the use of a 
track length estimator is to be preferred, as it would lead 
to a smaller variance. All analytical results have been 
validated by comparison with Monte Carlo simulations. 
For p < \ and transparent boundaries, the solutions of 



y m{(ti.)(xo)}, 



(32) 



for 771 > 1, where {-^^ and {-jy denote averaging Xq over 
the surface E (of V) or the volume V , respectively. Here 
d = 1 and m — 2, and it is easy to verify that 



{(t^)(xo)}, 



R 



{(4)(xo)}^ = 2(^i?^ + i 



R+-R^ 



R' 



(33) 



hence {(i^)(xo)}^ = {(t2,)(xo)}^ /2 {(ti.)(xo)}j.. Re- 
mark that we have an overall factor 1/2 with respect 
to the surface averages in [8l, because trajectories are 
here allowed starting from E in the outward direction, 
whereas in [5] they are not. In Fig. [2] we display the mean 



the rod problem ( 27 ) and ( 28 1 are given by combinations 
of exponential functions, rather than linear polynomials. 
In this case, the expressions for the moments are rather 
cumbersome and will not be reported here. Instead, we 
plot the moments as a function of the initial condition Xq, 
the domain size R and the scattering rate p. In Fig. |4] 
we display the mean collision number (ny)(xo) and the 
mean residence time (ty)(xo) for transparent boundaries 
andp = 0.5. The two surfaces, as a function of Xq and R, 
coincide, and this relation holds for any value of p. The 
second moments (ny)(xo) and (ty)(xo) are displayed in 
Figs. [5] (p = 0.9), [6] (p = 0.5), and[7](]5 == 0.1). In this 
case, it is not possible to establish a simple inequality 
between the two surfaces, independent of p. As the scat- 
tering rate varies, the surfaces change and there exist a 
value of p for which (ny)(xo) is smaller than (iy)(xo). 
This means that in presence of absorption the collision 
estimator may lead to a smaller variance. In Fig. [8] we 



collision number (ny)(xo) and the mean first-passage display the difference (ny)(xo) — (ty)(xo) when Xq = 0, 



time (iy)(xo) for leakage boundary conditions. The 
two surfaces, as a function of Xg and R, coincide, as 
expected from the considerations exposed above. This 



for various values of p: when R is large, (iy)(xo) becomes 
larger than (ny)(xo), and this behavior is enhanced for 
small values of p, i.e., large absorption rates. When R 



is large, the dependence on the angular variable lo gets 
progressively weaker, so that the integrals (201 and (14) 



coincide. Under this assumptions, calculations show that 
(ni.)(xo) ~ 1/(1 -p), and (n2,)(xo) ~ {l+p)/{l-p)\ 
independent of the initial condition. Then, from the 
equality of the Kac convolution integrals, the difference 
(ny)(xo) — (ty)(xo) for large R converges to the limit 
l/(p — 1), which implies a smaller variance for the colli- 
sion estimator. We have also verified that (iy)(xo) and 
(ty)(xo) satisfy the surface- and volume-averaged recur- 
sion property in [S], which generalizes Eq. (32) to resi- 
dence times. Again, all analytical results have been vali- 
dated by comparison with Monte Carlo simulations. 



IV. CONCLUSIONS 

Motivated by their relevance for stochastic transport 
phenomena as well as for Monte Carlo methods, in this 
paper we have examined the moments of collision num- 
ber nv and residence time ty of exponential flights in a 
volume V . We have presented a general approach that 
- based on the Kac moments formula - allows explicitly 
evaluating such quantities in terms of repeated convo- 
lutions of the particle equilibrium distribution. To ex- 
emplify the proposed formalism, we have in particular 
analyzed a Id system, the so-called rod problem, where 
closed form expressions can be found. We have therefore 
explicitly computed the moments of ny and ty for var- 



ious boundary conditions, focusing in particular on the 
first and second moment. Finally, the relevance of these 
findings in the context of Monte Carlo collision and track 
length estimators has been discussed. Results show that 
the averages of ny and ty coincide, whereas the second 
moments (hence the variances) depend on boundary and 
initial conditions, and on the scattering probability. Res- 
idence time has in general a smaller variance, but the op- 
posite is true when absorption dominates over scattering. 
By virtue of the increasing power of Monte Carlo meth- 
ods in solving realistic three-dimensional transport prob- 
lems, one might argue that such a simple system as the 
rod problem has a limited interest. On the contrary, 
we are persuaded that this analysis is useful, in that it 
allows focusing on the essential features of the physical 
system at hand. Indeed, on one hand it sheds light at 
the deep connections between sojourn times and collision 
number for exponential fiights, and on the other hand 
it gives some hints on the behavior of the intrinsic vari- 
ance of Monte Carlo collision and track length estimators. 
Extending the proposed approach to higher-dimensional 
and more complex systems is highly desirable, and inves- 
tigations to this aim are ongoing. 
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